Tunable disorder in a crystal of cold polar molecules 
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There is a growing interest in designing tunable many- 
body quantum systems, which can be used as physical 
simulators to access novel phenomena and to study the 
properties of analagous, but less controllable, many-body 
quantum systems [1-4]. For example, a major thrust of 
current research is focussed on realizing lattice-spin mod- 
els with ultracold molecules trapped on optical lattices 
[2, 3, 5-7]. This quantum simulation work relies on tun- 
ing interactions of ultracold molecules by external fields. 
While two-body interactions of ultracold molecules can 
be tuned by a variety of methods [8], it is much harder 
to achieve external field control over many-body interac- 
tions. In the present work, we demonstrate the possibility 
of controlling by an external field the dynamics of collec- 
tive excitations (excitons) of molecules on an optical lat- 
tice. We show that a suitably chosen two-species mixture 
of ultracold polar molecules loaded on an optical lattice 
forms a phononless crystal, where exciton-impurity inter- 
actions can be controlled by applying an external electric 
field. This can be used for the controlled creation of 
many-body entangled states of ultracold molecules and 
the time-domain quantum simulation of disorder-induced 
localization and derealization of quantum particles. 

The present work builds on three major recent ad- 
vances in atomic, molecular and optical physics: the cre- 
ation of dense ensembles of ultracold polar molecules in 
the ro- vibrational ground state by photoassociation of al- 
kali metal atoms [9, 10]; the development of techniques 
for trapping ultracold atoms and molecules by optical 
lattices [11]; and the development of techniques for de- 
tecting and manipulating single atoms in optical lattices 
[12-14]. It has thus become technologically possible to 
create ordered ensembles of ultracold polar molecules in 
the ro- vibrational ground state trapped by an optical lat- 
tice with a lattice separation of about 400 nm [15]. Using 
the method of Ref. [12], molecules in specific lattice sites 
can be evaporated. By overlapping two optical lattices 
of trapped molecules or using the method of Ref. [16], 
it should be possible to create a mixture of molecules 
arranged in an arbitrary array with any spatial dimen- 
sionality. 

Figure 1 illustrates the system considered here. We 
consider an optical lattice of 1 Y, polar molecules, such 
as LiCs, with one molecule in the rovibrational ground 
state per lattice site. We assume that the tunneling of 
molecules to different lattice sites is entirely suppressed, 
which can be easily achieved by applying laser fields of 
sufficient power [11]. The lattice sites are separated by 




FIG. 1: Schematic illustration of a crystal with tunable impu- 
rities. Molecules are confined by the periodic potential of an 
optical lattice and coupled by the dipole-dipole interaction, 
which allows the exchange of a rotational excitation between 
molecules in different lattice sites. An impurity is a molecule 
with a different dipole moment and different rotational energy 
splitting. For certain combinations of molecules, the energy 
defect Vb can be tuned by an external electric field from zero 
to a large magnitude, thereby modifying the energy transfer 
in the crystal. 



400 nm so the interaction between the molecules is en- 
tirely determined by the long-range dipole-dipole inter- 
action potential. We consider the lowest excited states 
of the molecular crystal. In the absence of the dipole - 
dipole interactions, these would be the rotational exci- 
tations of the individual molecules. However, the inter- 
actions between the molecules lead to collective modes 
known as Frenkel excitons [17]. Each exciton state is a 
plane wave characterized by a wavevector. In the ab- 
sence of an electric field, the rotational ground state of a 
diatomic molecule is characterized by the rotational 
angular momentum N — 0. The first excited rotational 
state N — 1 is triply degenerate. The N = — > N = 1 
excitation of molecules therefore leads to three excitonic 
states characterized by the dispersion curves shown in 
Figure 2. In the presence of an electric field, the rota- 
tional state N = 1 splits into two Stark energy levels 
corresponding to the projections M = and M = ±1 
of N on the electric field axis. This lifts the degeneracy 
of the exciton states as shown in Figure 2a. When the 
electric field is perpendicular to the molecular array, ex- 
citons 7 and f3 have negative effective mass, i.e., their 
energy decreases with increasing wavevector, while the 
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FIG. 2: Stark effect on rotational excitons. Exciton energy 
E(k) calculated for three lowest excitations in a ID array of 
10 5 LiCs molecules separated by 400 nm as functions of the 
exciton wavenumber k. Exciton states a, (3 and 7 centered 
at the isolated molecule transition energies AE eg are shown 
in the presence of an electric field perpendicular to the array 
(panel a) and parallel to the array (panel b). The electric 
field magnitudes are 2 kV/cm (dotted lines), 3.2 kV/cm (solid 
lines), 5 kV/cm (dashed lines). At these electric fields, exciton 
7 is separated from excitons a and /? by the energy much 
larger than the exciton band widths. The figure shows that the 
magnitude and the sign of the effective mass of the rotational 
excitons can be controlled by varying the magnitude and the 
direction of the electric field. 

effective mass of exciton a is positive. The effective mass 
can be changed by varying the strength and the direc- 
tion of the external electric field (see Figure 2b) . We also 
note that, unlike electronic excitons in molecular crystals, 
rotational excitons have a very long radiative lifetime ex- 
ceeding seconds. In the following, we focus on the exciton 
state labeled by 7 in Figure 2a. 

If a small number of LiCs molecules are replaced with 
molecules of a different kind, such as LiRb, the trans- 
lational symmetry of the lattice is disturbed, as in a 
solid crystal with impurities (see Figure 1). Let AE eg 
be the rotational excitation energy in host molecules, 
AE eg + Vq the energy of the same transition in impu- 
rities, and AJ n ^ m the difference between impurity-host 
and host-host coupling constants for molecules at sites n 
and to. Consider first the lattice with a single impurity 
at site n = 0, which is described by the Hamiltonian [18] 

H = H + VoBlB + A Jn,o (bI Bo + B n B^j , (1) 

where i?o is the free-exciton Hamiltonian and Bq and 
Bo are the exciton creation and annihilation operators 
defined in the Methods section below. The exciton- 
impurity interaction can be described as a <5-function po- 
tential with strength Vq, and a perturbation due to the 
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FIG. 3: Excitation energies of non-inetracting molecules in an 
electric field. Panel (a) shows the excitation energies AE eg 
for transitions \N = 0,M = 0) \N = 1,M) with M = 
(upper curve) and M — ±1 (lower curve) vs electric field for 
three polar molecules. Panel (b) shows an expanded view of 
the encircled area in panel (a). Panel (c) shows the exciton- 
impurity 2D scattering cross sections for exciton 7 in Figure 
(2a) with |k|a = 4 x 10" 5 (solid line), |k|a = 4 x 10" 3 (dashed 
line), and |k|a = 4 x 10 -2 (dotted line). The calculations are 
for an array of LiCs molecules with one LiRb impurity. The 
lattice constant is a = 400 nm and £0 — 3228.663 V/cm. 

difference in the dipole moments of host and impurity 
molecules. 

We propose to modify V by shifting the rotational lev- 
els of host and impurity molecules simultaneously using 
a static electric field. Figure 3a shows the excitation en- 
ergies AE eg for the rotational transition N = — > 1 of 
CsF, LiCs and LiRb molecules as functions of the elec- 
tric field and demonstrates that the excitation energies 
for these combinations of molecules become the same at 
certain values of the electric field. The interaction of 
polar molecules with the electric field allows for the pos- 
sibility to explore the dependence of the exciton-impurity 
scattering cross section not only on the exciton wavevec- 
tor k, but also on the scattering strength V , which is 
not possible in conventional solids. Figure 3c shows the 
cross section for elastic scattering of an exciton by a sin- 
gle impurity in a 2D crystal as a function of the electric 
field, demonstrating resonant enhancement of the scat- 
tering cross section occurring for values of Vq that sup- 
port a shallow bound state (see Methods). These reso- 
nances are analogous to Feshbach resonances in atomic 
collisions, commonly used to tune the scattering prop- 
erties of ultracold atoms [19]. Figure 3c shows that the 
scattering properties of excitons in the system proposed 
here can be tuned by varying the external electric field. 

Consider now a lattice with multiple impurities. It is 
well known that quantum particles in the presence of a 
random distribution of scattering centers undergo coher- 
ent localization [21]. In the system proposed here, the 



3 



localization of excitons can be tuned by an electric field. 
As an example, we consider localization of excitons in a 
ID array of LiCs separated by 400 nm, with a random 
homogeneous distribution of LiRb impurities (see Figure 
4). The Hamiltonian describing an exciton in the pres- 
ence of Ni substitutional impurities at positions \ n can 
be written as H = H + W + V, with the corresponding 
matrix elements in the basis of free-exciton states |^k) 
given by 



(#o) q ,k = E(k)6 k . 
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where E(k) is the energy of the exciton, k and q de- 
note exciton wavevectors and AJ(a) = AJ n n _ 1 . If the 
impurities are distributed randomly, the terms V and W 
correspond to diagonal and off-diagonal disorder, respec- 
tively. 
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FIG. 4: Probability density |^(x)| 2 describing an exciton 
near the top of the energy spectrum for a ID array of 1000 
LiCs molecules with 10% of homogeneously and randomly dis- 
tributed LiRb impurities. Panels correspond to different val- 
ues of V : (a) Vb = 0, (b) Vo/h = 22 kHz, and (c) Vo/h = 100 
kHz. Tuning Vb by an electric field can thus be used to realize 
different regimes of exciton dynamics from free propagation 
to strong localization in different spatial arrangements. The 
difference of the dipole moments of LiCs and LiRb molecules 
leads to the value A J = —6.89 kHz. 

When an external electric field is such that Vq = 0, 
exciton-impurity scattering occurs only due to the 
difference in dipole moments between host and impurity 



molecules. Figure 4 shows the probability density of a 
particular eigenstate of Hamiltonian (2) near the top 
of the energy spectrum, for different values of Vq. Due 
to the negative effective mass of free-exciton states, 
high energy eigenstates are dominated by free-exciton 
states with k s» 0. These eigenstates are localized (Fig. 
4a). Delocalization of these states can be achieved 
by applying an electric field so that Vq m — 4AJ(a) 
(see Fig. 4b and Eqs. (2)). In this case, for a given 
k, the matrix elements (V)q t u and (W)q } u cancel for 
q ps k, which suppresses the coupling between the 
corresponding free-exciton states. Localized eigenstates 
with different energies become delocalized at different 
values of Vq. The wavepackets are localized for values 
of Vo that do not balance the effect of A J (a) (Fig. 4c). 
Figure (4b) shows an example of delocalization in a ID 
disordered system due to correlations between diagonal 
and off-diagonal disorder [20]. 

Applications 

The system proposed here offers three unique features: 
(i) long-lived excitons that are stable against sponta- 
neous decay and whose effective mass can be controlled 
by an external electric field; (ii) dynamically tunable im- 
purities; (iii) possibility to arrange impurities and host 
molecules in various configurations and dimensions. This 
opens up new possibilities for quantum simulation of fun- 
damental physical phenomena: 

Exciton - impurity scattering leads to localization of 
excitons in crystalline solids [21]. The exciton - im- 
purity interactions in the system proposed here can be 
tuned from resonantly enhanced (Figure 3c) to entirely 
suppressed (Figure 4b) scattering, which can be used to 
study the dynamics of exciton localization, including the 
timescales for the formation of localized states and their 
dependence on exciton bandwidth and exciton - impurity 
interaction strength. 

The negative effective mass of the excitons may be used 
for studies of negative refraction of microwave fields. At 
present, the materials exhibiting negative refraction are 
artificial metamaterials [22], in which the possibility to 
use a macroscopic (spatially averaged) description of elec- 
tromagnetic fields is questionable [23] . We also note that 
the localization of quasi-particles with negative effective 
mass, which has not yet been observed, may exhibit dif- 
ferent mechanisms of quantum interference due to mul- 
tiple scattering events. Measuring the dynamics of exci- 
tons with tunable effective mass may provide new insights 
into the mechanism of Anderson localization [21]. 

Controlled spatial distributions of impurities and 
molecular crystals with specific arrangements of crystal 
particles in one, two or three dimensions may be used 
to study the effects of dimensionality and finite size on 
energy transfer in mesoscopic materials [24]. In addi- 
tion, the system proposed here is ideally suited for the 
study of long- and short-range correlations in the disor- 
dered potentials. The presence of short- and long-range 
correlations in the disordered potentials may result in 
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the appearance of a discrete [25, 26] or even continu- 
ous [27, 28] set of dclocalized states in low-dimensional 
crystals. Measurements of exciton localization in a crys- 
tal with tunable impurities can be used for time-domain 
quantum simulation of disorder-induced localization and 
derealization of quantum particles. 

The localized states displayed in Figure 4 are many- 
body entangled states of the molecules in an optical lat- 
tice. The possibility to tunc exciton - impurity interac- 
tions can thus be exploited for the controlled creation of 
many-body entangled states of ultracold molecules, nec- 
essary for the experimental realization of quantum com- 
putation with molecular ensembles [29]. 

Finally, a mixture of ultracold molecules with impu- 
rities forming a sublattice may be used to study the 
formation of wavevector space crystals of excitons. The 
eigenstates of such a two-species lattice correspond to a 
discrete distribution in k-spacc. Tuning the impurities 
by a sinusoidally varying electric field could then be 
used to induce resonant transitions between different 
fc-states. In particular, it could be interesting to couple 
reversibly the low-fc optically active, but slow (small 
group velocity) states with optically inactive, but fast 
(larger group velocity) states from the middle of the 
Brillouin zone. This would lead to the creation of 
excitonic wave packets, which can be slowed down or 
accelerated by an external electric field. 

Methods 

Dispersion curves 

The Hamiltonian describing the optical lattice with iden- 
tical molecules in the presence of a static electric field is 



H = 



iV mo i 

E 

71=1 



(B e Nl - d„ 



• E 



n—1 rriy^n 



V dd (r n - r TO ), 



(3) 

where r„ is the position of the n-th lattice site, N n is the 
rotational angular momentum of the molecule, B e is the 
rotational constant, d„ is the electric dipole operator, E 
is a static electric field, V dd is the dipole-dipole interac- 
tion between molecules in different lattice sites, and N mo \ 
is the total number of molecules. Hamiltonian (3) can be 
written as [17] 



1 



H - AE eg ^ B l,M B n,M + 2 E J n,m B n,M B m,M 



1 



+ 2 2^ E J n,m B l,M B m,M'- (4) 

The exciton creation operators arc defined by 
B n,M\9)n = \e,M) n , where \g) and |e, M) denote 
ground and excited field-dressed rotational states, 
which have a well-defined projection M of the rota- 
tional angular momentum TV along the direction of E. 
The excitation transfer between molecules at lattice 
sites r„ and r m is described by the coupling constant 



have neglected the gas - condensed matter shift in the 
diagonal terms since it is much smaller than A.E eg . 

If E is perpendicular to the intermolecular axis R = 
fn — r m, then the dipole-dipole interaction between 
molecules n and m couples two-molecule states with 
A(M„ + M m ) = 0, ±2, and if E is parallel to R, the 
selection rule is A(M„ + M m ) = 0. 



By using B\ 



t £ 



* kr n where k is 



the exciton wavevector, Hamiltonian (4) can be rewritten 

as 

H = E [^^eg^M,M' + L M ,M'(k)] B^ M B^ M i, 

k M.M' 



T M,M' ik-r„ 
Z^n J n,0 e 



(5) 

Hamiltonian (5) 



with L M ,M'{k) 

is diagonalized by the unitary transformation & 



k.M 



y ] ctM^A^ u [17] . The resulting Hamiltonian is H 



k.^i 

Ek^-^W^k.^k^, where fi = a,f3,j labels the dis- 
persion relations associated with the eigenstates. 

In the case where E is parallel to the molecular array, 
the eigenstates of Hamiltonian (5) are: |ak) = |*k,M=i), 
I Ac) = l*k,M=-i), and | 7k ) = |* k , M=0 ). If the electric 
field E is perpendicular to the molecular array, the eigen- 
states are: |a k ) = ^ [|*k,M=i) - |*k,M=-i>], \Pk) = 
^ [|*k,M=i> + |*k,M=-i>], and | 7k ) = |* k ,M=o)- 
The dispersion curves ^(k) of the three exciton 
branches are shown in Fig. 2 for a ID array of LiCs 
molecules (d = 5.529 Debye, 2B e = 11.7998 GHz 
[30]) with lattice constant a = 400 nm in the presence 
of an electric field parallel and perpendicular to the 
array. Analogous dispersion curves can be obtained 
in 2D and 3D and for other directions of the electric field. 

Scattering resonances of excitons 

Figure 3c shows the 2D exciton-impurity elastic cross sec- 
tion as a function of the electric field, obtained in the 
effective mass approximation E{k) = E(0) — T^, for neg- 
ative effective mass. = fi 2 |k| 2 /2|m*| is the kinetic en- 
ergy of the exciton. The impurities are LiRb molecules 
with d = 4.165 Debye, 2B e = 13.2268 GHz [30]. If E b 
is the energy of the bound state produced by Vo, and A 
is the exciton bandwidth, we write the scattering cross 
sections for |k| <C n/a, in 3D and 2D, as 



c 3 D(fc, V ) 



27rft 2 
T k + E, 



<72d(&, Vq) 



4lT 2 /k 



n 2 + In 2 



; (A-T fc 



Tk(A — E b ') 



(6a) 



(6b) 



M.M' 



(e,M\ n (g\ m V dd \g) n \e',M') m . In Eq. (4), we 



where E^ D) = (2/tt - V Q /A) 2 A and E ( ^ D) = 
A/ [exp(4A/7rVo) — 1]. In 3D, the bound state exists 
only if V > 2A/ir. In 2D, any positive V < 0.2 A 
supports an exponentially shallow bound state. Reso- 
nant enhancement of the scattering cross section occurs 
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for values of Vq that support a shallow bound state 
(Eb <~ T — > 0). In ID, shallow bound states appear 
only for vanishingly small Vq and resonances may not 
be observable. As a result of the negative effective mass 
of the exciton, the bound state that leads to resonance 
is produced by a repulsive potential. Equations (6) are 
derived in the approximation AJ„ = 0. Including 
A J n fi in the calculation leads to a shift of the positions 



of the resonance [18] and the resonant enhancement of 
the scattering cross section at a slightly different value 
oiV Q . 
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